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EQUATIONS  OF  MOTION  FOR  NONAXISYMMETRIC 
VIBRATIONS  OF  PROLATE  SPHEROIDAL  SHELLS 


1.  INTRODUCTION 

Prolate  spheroidal  shells  (PSS)  are  shells  of  revolution  that  can  span  a  wide  range  of 
shapes,  from  spherical  to  needle-like.  They  are  essentially  ellipsoids  of  revolution  about  their 
major  axis  and  may  have  a  constant  thickness  or  a  variable  thickness  defined  by  two  confocal 
ellipses.  The  ellipse  is  characterized  by  an  interfocal  length  d  and  eccentricity  1/a,  where  a  is  the 
radial  coordinate  of  the  shell's  mid-surface  in  the  prolate  spheroidal  coordinate  system.  In  the 
limits  d  ->  0,  a  -»  and  (da/2)  ->  R,  the  ellipse  approaches  a  circle  whose  radius  is  R  and  the 
shell  then  becomes  spherical. 

The  theory  of  thin  elastic  shells  is  based  on  the  approximation  that  the  shell  thickness  h  is 
much  smaller  than  the  typical  dimension  of  the  shell  R,  i.e.,  h/R  «  1.  Typically,  the  ratio  of 
h/R  is  taken  to  be  less  than  0.05.  Unlike  thin  elastic  plates,  the  in-plane  and  out-of-plane 
displacements  of  thin  elastic  shells  are  coupled.  Furthermore,  unlike  elastic  plates,  elastic  shells 
can  deform  in  an  extensional  mode  (membrane)  or  in  a  coupled  extensional/flexure  mode 
(membrane  and  bending). 

In  the  extensional  theory  of  thin  shells  of  revolution,  there  are  three  independent 
displacements:  w  is  the  out-of-plane  displacement,  u  is  an  in-plane  displacement  (extension)  in  the 
plane  of  the  generator  surface,  and  v  is  a  displacement  perpendicular  to  the  generator  surface,  i.e., 
torsional.  The  displacements  u  and  w  are  coupled  in  longitudinal  and  flexural  deformation,  while 
v  is  uncoupled  and  is  purely  torsional  deformation.  Therefore,  the  strain  energy  density  is 
exclusively  dependent  on  the  extensional  strains  of  the  shell's  mid-surface  and  is  linearly 
dependent  on  the  thickness  h/R.  Thus,  for  the  extensional  theory  of  shells,  there  are  three 
branches  in  the  frequency-wavenumber  spectrum  corresponding  to  the  three  displacements  u,  v, 
and  w.  The  lowest  branch  is  an  acoustic  branch  representing  flexural  vibrations  of  the  shell  where 
|w/u|  >  1.  This  acoustic  branch  has  a  constant  limit  in  the  frequency  as  the  wavenumber  k  «>. 
This  means  that  the  group  and  phase  velocities  for  this  branch  approach  zero  as  k  °°,  rather  than 
the  sound  speed  of  the  shell  material.  Furthermore,  the  flexural  resonance  frequencies  are  not 
distinct  for  k  since  the  frequency  spacing  approaches  zero  as  the  mode  number  k  becomes 
very  large.  The  theory  predicts  the  low-ordered  mode  resonance  frequencies  accurately,  but 
becomes  inaccurate  after  the  first  few  modes.  The  second  acoustic  branch  represents  a  decoupled 
(pure)  torsional  displacement  v,  which  has  group  and  phase  velocities  correctly  approaching  the 
shear  sound  speed  of  the  shell  material,  i.e.,  cs  =  G  /  ps  as  k  ->  where  G  is  the  shear  modulus. 
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and  ps  is  the  density.  The  third  branch  is  an  optical  branch  that  represents  |w/u|  <  1,  i.e.,  primarily 
in-plane  deformation  of  the  shell,  representing  longitudinal  vibrations.  As  k  ->  °°  this  branch  has 
phase  and  group  velocities  equal  to  the  sound  speed  in  plate  materials  (i.e.,  c2  =E/(ps(l  -  v2)), 
where  E  is  the  Young's  modulus  and  v  is  the  Poisson's  ratio),  which  are  close  to  the  exact  values. 

In  the  coupled  extensional-bending  theory  of  thin  shells,  h/R  «  1  is  again  assumed. 
However,  the  change  of  the  curvatures  of  the  shell’s  mid-surface  and  (32  are  now  included  in 
the  strain  energy  density  of  the  shell.  These  changes  of  curvature  (bending)  are  expressed  in 
terms  of  the  derivatives  of  u,  v,  and  w,  so  that  the  strain  energy  density  is  again  dependent  on 
only  these  three  displacements.  The  strain  energy  has  membrane  components  that  are  linearly 
dependent  on  the  thickness  (h/R),  and  additional  bending  components  have  cubic  dependence  on 
the  thickness,  i.e.,  (h3/R3).  Thus,  there  are  still  three  branches  in  the  frequency-wavenumber 
spectrum  corresponding  to  the  three  independent  variables  u,  v,  and  w.  The  lower  acoustic 
branch  representing  the  flexural  vibration  is  dependent  on  terms  of  (h/R)  and  (h3/R3).  However, 
the  group  and  phase  velocities  become  infinite  as  k  ->  i.e.,  to  <=  k^,  similar  to  the  Bemoulli- 

Euler  bending  theory  of  elastic  plates.  Thus,  the  flexural  vibration  is  accurate  in  the  low 
frequency  range  but  is  still  inaccurate  in  the  high-frequency  range.  The  second  acoustic  branch 
is  the  extensional  torsional  mode,  which  is  not  influenced  by  the  addition  of  bending.  The  optical 
branch  representing  longitudinal  vibrations  is  slightly  influenced  by  the  bending  terms;  and  thus, 
the  group  and  phase  velocities  are  good  approximations  to  the  exact  values  in  the  limit  of  k  ->  °°. 

To  improve  the  accuracy  of  the  flexural  vibration  branch  for  high  frequencies,  the  thin  shell 
theory  is  improved  by  the  addition  of  shear  deformations  and  rotatory  inertias.  This  means  that 
the  change  of  curvatures  (5j  and  P2  are  no  longer  derivatives  of  u,  v,  and  w,  but  are  independent 
variables.  The  new  thin-shell  theory  now  has  five  independent  variables,  (u,  v,  w,  P2,  and  (32); 
the  frequency-wavenumber  spectrum  has  five  branches:  two  acoustic  branches  (one  torsional  and 
one  flexural)  and  three  optical  branches  (one  longitudinal  and  two  thickness  shear  modes).  The 
flexural  branch  now  has  the  correct  group  and  phase  velocities. 

This  report  first  surveys  the  literature  covering  prolate  spheroids  and  PSS  for  free  and 
forced  vibrations  of  unloaded  shells,  fluid-filled  shells,  free  and  forced  vibrations  of  submerged 
shells,  and  acoustic  scattering  from  submerged  shells.  It  then  presents  the  development  of  the 
five  coupled  partial  differential  equations  for  thin  elastic  shells,  including  extensional,  bending, 
shear  deformation  and  rotatory  inertias.  These  equations  are  developed  from  the  Lagrangian  of 
the  strain  and  kinetic  energy  densities  of  the  shell,  using  Hamilton's  variational  principle.  This 
approach  leads  to  an  accurate,  self-consistent  coupled  set  of  five  partial  differential  equations  of 
the  correct  expansion  in  h/R  valid  to  0(h3/R3).  This  approach  also  guarantees  a  self-adjoint 
positive-definite  system  of  equations  that  leads  to  real  and  positive  eigenvalues  and  real 
eigenfunctions. 
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2.  LITERATURE  REVIEW 


2.1  FREE  AND  FORCED  VIBRATIONS  OF  UNLOADED  PROLATE 
SPHEROIDAL  SHELLS  (PSS) 

DiMaggio  and  Silbiger1  were  the  first  to  examine  the  free  extensional  axisymmetric 
torsional  vibrations  of  a  PSS  of  variable  thickness.  The  equations  of  motion  were  based  on  the 
Lagrangian  of  the  shell  using  Hamilton's  variational  principle.  The  exact  solution  was  found  in 
terms  of  the  angular  prolate  spheroidal  wave  function  SonCn).  Silbiger  and  DiMaggio2  developed 
the  equations  for  axisymmetric  extensional  flexural  vibration  of  PSS  of  a  variable  thickness. 
Using  the  Rayleigh-Ritz  method,  they  obtained  the  resonance  frequencies  and  mode  shapes,  but 
the  solution  was  limited  to  low  eccentricity  of  the  shell  (  =  0.7).  Shiraishi  and  DiMaggio3 
obtained  the  solution  by  using  perturbation  techniques.  The  perturbation  parameter  was  the 
eccentricity  (1/a)  and  their  perturbation  series  was  extended  to  0(l/a4).  It  should  be  noted  that 
the  zero  order  term  is  the  spherical  shell  resonances  and  mode  shapes.  Again  these  frequencies 
were  valid  for  low  eccentricity  ( =  0.7).  Nemergut  and  Brand4  developed  the  extensional  flexural 
vibration  of  a  constant-thickness  shell  using  a  numerical  integration  scheme.  DiMaggio  and 
Rand5  obtained  the  solution  of  variable-thickness  flexural  extensional  vibration  by  finite 
difference  techniques.  Rand6  obtained  the  exact  solution  for  torsional  vibration  of  solid  prolate 
spheroids  and  prolate  spheroidal  shells  using  the  exact  theory  of  elasticity. 

Burroughs  and  Magrab7  were  the  first  to  include  bending  and  shear  deformation  in  the  shell 
equations  using  a  variational  approach  in  terms  of  shear  and  moment  resultants.  They  developed 
five  coupled  equations  for  nonaxi symmetric  motion  for  a  shell  with  constant  thickness. 

However,  they  obtained  solutions  for  only  the  axisymmetric  vibration  using  Galerkin's  method, 
which  means  that  there  are  three  branches  for  the  axisymmetric  case.  They  showed  that  for 
axisymmetric  vibrations,  the  resonance  frequencies  for  the  flexural  modes  increase  with  (h/R), 
while  the  longitudinal  modes  are  not  influenced  by  the  addition  of  bending  or  shear 
deformations.  Of  course,  there  is  a  new  branch  representing  thickness-shear  mode  that  was 
shown  to  be  inversely  proportional  to  (h/R).  Yahner  and  Burroughs8  obtained  the  axisymmetric 
vibration  of  an  open  shell  with  a  variable  thickness.  They  exhibited  the  influence  of  the  opening 
of  the  shell  on  the  three  axisymmetric  modes,  flexural,  longitudinal,  and  thickness-shear.  Chen 
and  Ginsberg9  revisited  the  problem  by  focusing  on  the  loci  of  the  eigenvalues  of  three  PSS  for 
the  axisymmetric  vibration,  but  they  included  only  extensional  and  bending  terms  and  used  the 
Galerklin  method  in  their  solution. 
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2.2  VIBRATION  OF  FLUID-LOADED  PSS 


The  first  attempt  at  examining  the  fluid  loading  effects  on  the  vibration  of  submerged  PSS 
was  made  by  Hayek  and  DiMaggio. 10,11  The  complex  natural  frequencies  of  submerged  PSS  using 
only  extensional  strain  energy  were  found  for  shells  of  variable  thickness  using  perturbation 
techniques  on  the  strain  energy  density  and  the  surface  acoustic  pressure.  They  showed  that  the 
two  pairs  of  real  frequencies  (±o>i  2)  °f  the  flexural  and  longitudinal  branches  of  in- vacuo  shells 
have  increased  to  seven  frequencies:  three  pairs  of  complex  frequencies  (±{012,3  +  i£j  2  3)  and  one 
purely  imaginary  frequency  (i^).  The  two  pairs  with  small  imaginary  components  (±(Dj  2+  i^i  2) 
of  the  submerged  shells  correspond  to.  the  two  real  pairs  of  the  unloaded  shells.  The  real  parts 
(±g>i>2)  decreased  from  the  in-vacuo  shells  due  to  the  additional  fluid  mass  loading  of  the  acoustic 
medium,  and  the  imaginary  parts  represent  the  energy  lost  through  acoustic  radiation  (acoustic 
resistance)  from  the  shell.  The  remaining  pair  ( ±co3  +  i£3)  has  a  large  imaginary  part,  caused  by 
high  acoustic  resistance.  The  purely  imaginary  root  implies  pure  damping.  Rand  and  DiMaggio12 
examined  the  resonance  frequencies  of  extensional  vibration  of  fluid-filled  PSS.  The  presence  of 
a  fluid  inside  the  shell  increases  the  number  of  resonance  frequencies  because  of  the  three- 
dimensionality  of  the  cavity  resonances  of  the  interior  of  the  PSS.  Yen  and  DiMaggio13  obtained 
the  solution  for  the  forced  vibration  of  submerged  PSS  due  to  axisymmetric  time-harmonic  surface 
forces.  They  used  numerical  integration  techniques  after  transforming  the  infinite  exterior 
acoustic  medium  domain  to  a  finite  domain.  Bedrosian  and  DiMaggio14  obtained  the  transient 
solution  of  an  axisymmetric  extensional  vibration  of  a  submerged  PSS  undergoing  sudden 
application  of  a  uniform  surface  source  at  t  =  0. 

Berger15  was  the  first  to  include  bending  effects  in  the  problem  of  forced  vibration  of  a 
submerged  PSS  of  constant  thickness.  He  used  a  finite  difference  approach  for  the  shell 
equations  and  mapped  the  infinite  region  of  the  exterior  acoustic  medium  into  a  finite  region. 
Berger  presented  the  response  of  the  shell  attributed  to  a  normal  transient  loading.  Lee  and 
DiMaggio16  obtained  the  response  of  a  fluid-filled  PSS  as  a  model  of  a  human  head.  They 
included  bending  theory  in  their  model  and  obtained  the  response  of  the  PSS  due  to  harmonic 
surface  sources.  Ross  and  Johns17  obtained  the  solution  for  the  axisymmetric  vibration  of  a  free 
and  submerged  (both  sides)  hemi-PSS  and  conducted  experiments  to  support  their  analytic 
predictions.  Prikhod'ko18  obtained  the  solution  for  a  slender  submerged  PSS,  but  the  solution  is 
valid  for  sufficiently  slender  shells  in  order  to  neglect  the  fluid  loading  on  the  shell  ends. 
Pauwelussen19  used  a  finite  element  method  (FEM)  to  solve  the  problem  of  axisymmetric 
response  of  submerged  PSS  due  to  a  shock  load.  Chen  and  Ginsberg20  obtained  the  acoustic 
radiation  for  PSS  including  bending  effects  due  to  axisymmetric  mechanical  surface  forces. 
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2.3  ACOUSTIC  SCATTERING  FROM  PSS 


Acoustic  scattering  from  PSS  has  not  been  studied  as  well  as  the  topics  of  vibration  and 
acoustic  radiation  from  PSS.  Chertock,  Hirsh,  and  Ogilvie21  obtained  the  flexural  response  of 
PSS  to  incident  underwater  explosion  and  tested  the  theory  with  experiments  on  a  paralellopiped 
box-shell.  Silbiger22  studied  scattering  from  PSS  due  to  an  incident  time-harmonic  planewave  by 
approximating  the  response  from  one  resonant  mode.  Jones-Oliveira23'24  obtained  the  transient 
scattering  from  an  axially  incident  plane  wave*  i.e.,  axisymmetric  response,  from  a  PSS  that 
includes  extensional  and  bending  effects  using  expansions  in  terms  of  Legendre  polynomials. 

She  also  obtained  the  resonances  of  two  different  PSS  for  the  first  eight  flexural  and  longitudinal 
modes. 
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3.  PROLATE  SPHEROIDAL  COORDINATE  SYSTEM 


In  this  report,  the  prolate  spheroidal  coordinate  systems  of  Flammer25  and  Hanish,  Baier, 
Van  Buren,  and  King26  have  been  employed.  Basically,  these  coordinates  represent  ellipsoidal 
hyperboloidal  surfaces  (see  figure  3.1).  The  coordinate  system  transformations  from  cartesian 
coordinates  are  given  by 


x  =  —  [(1  -  T]2 )  (£2  - 1)]1/2  cos  0  -1  —  rj  —  1, 


y  =  ^[(1  -  vi2)(¥  ~  1)]1/2  sin*  1  <  %  < 


d  . 

Z  =  — TIC 
2 


0  <  *  <  2n, 


where  d  represents  the  interfocal  distance  of  the  generating  ellipse. 

In  terms  of  cartesian  coordinates,  the  prolate  spheroidal  coordinates  are 


2,2  2 

X  +  V  z 

— 2 — - + - 2 —  =  *  (ellipsoid  surface  for  constant  |), 

I)  «•->  (I)  5’ 


2  2  2 
y*r  +  y  z 

—7 - - - -j - =  -1  (hyperboloid  surface  for  constant  rj), 

(|)2(  l-Tl2)  (~)V 


*  =  arctan  (y/x)  (half-plane  surface  for  constant  *). 


The  ratio  of  minor  to  major  axes  of  a  prolate  spheroid  of  constant  £,  is 


_SL l 


\>  L 
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Figure  3.1.  Prolate  Spheroidal  Coordinate  System 


and  the  eccentricity  of  the  prolate  spheroid  is 


1 


Thus,  the  prolate  spheroid  represents  many  ellipsoidal-shaped  objects,  e.g., 
£  ->  1,  a  =  0,  a  straight  line  of  length  d, 

£  °°,  cx  =  1,  a  sphere  of  constant  radius  =  (d£)/2,  d  ->  0,  |  «>, 

T)  =  ±1,  semi-infinite  lines  along  the  z-axis. 


rj  =  constant,  which  represents  a  hyperboloid  that  makes  an  asymptotic  angle  0  =  arccos  (r\) 
with  the  z-axis,  and 

r)  =  0,  the  x-y  plane. 

Note  that  in  the  limit  t,  °°,  d  ->  0,  the  prolate  spheroidal  coordinates  reduce  to  spherical 
coordinates,  i.e.,  d  t/2  r,  r)  ->  cos  0,  as  d  0  and  £  -> 

The  length  and  area  elements  are 


ds2  =  dx2  +  dy2  +  dz2  =  h^  dr)2  +  h|  d£2  +  h|  d<|>2’ 


where 


hn  =  — 
2 


(t  2  Jl\ 


¥~r\ 

i 

1-Tl 


1/2 


’  Vi 


(tl 


r-i 


1/2 


,  and  h^=^[(l-riz)(r-l)]1/2. 


(3.5) 


An  area  element  on  a  surface  t,  =  constant  is  defined  as 


dA  =  h^h^dri d<j)  =  (^-J  ^  -ti2)(|2  -1) dr) d<(> , 


(3.6) 


and  the  gradient  is  defined  by 


n,lf=  IjETaw  ll-ri2  3\|/  1  IjV 

V  (d/2)^2-n2  dt,  (d/2)^2_n2  0T1  (d/ 2) -y (^2  — 1)(1  —  n2)  3<i>  ' 


(3.7) 


The  Helmholtz  equation  in  prolate  spheroidal  coordinates  is 


(V2  +k2)\|/  =  0 


j-a-'n2)^  jr(£2  -1)^+ r,2  ~  ~ — 27^t +c2(^2-,n2) 
01)  3r|  dt,  dt,  (£2-l)(l-if)a<r 


V  =  0,  (3.8) 


where  the  nondimensional  wavenumber  c  is  given  by 


c=  -kd. 
2 


(3.9) 
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The  Helmholtz  equation  is  separable  into  prolate  spheroidal  wave  functions,  i.e., 

~cos(m<|>)~ 

sin(m<|>) 


Vmn  ~  ^mn(c>1l)^mn(C’^) 


(3.10) 


where  the  angular  and  radial  wave  functions  Smn(c,ri)  and  Rmn(c£)  satisfy  the  following  second- 
order  differential  equations: 


_d_ 

dri 

_d_ 

d£ 


(1  ^  Smn  (c,  T|) 


«2-l)-^Rmn(c,Tl) 


y  „2„2  m 

Amn  “c  1}  ~ 


i — ry^ 


d£ 


Kn~C2Z>2  + 


m 


sz-i 


Rmn(^’^l)  —  0, 


(3.11) 


(3.12) 


where  is  the  separation  constant. 


Note  that  if  t,  d  ->  0,  c  0,  £  d/2  r,  c£  kr,  and  T|  ->  cos  0,  then  these  equations 
reduce  to  the  wave  equation  in  spherical  coordinates.  Thus,  Xmn  ->£(£  +  1),  where  £  =  m,  m  + 
1,  m  +  2, ...  and 

Smn(c>  Tl)  P n  0l)  ar)d  Q  “(tl)  (associated  Legendre  functions), 

Rmn(c,  4)  *  jnCkr)  and  yn(kr)  or  h  ^(kr)  and  h  (kr)  (spherical  Bessel  functions). 

It  should  be  noted  that  the  angular  prolate  spheroidal  wave  functions  are  dependent  on 
wavenumber  c  unlike  the  angular  spherical  wave  functions.  Both  of  the  prolate  spheroidal  wave 
functions  can  be  expanded  in  terms  of  spherical  wave  functions. 
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4.  DEFORMATIONS  AND  KINEMATICAL  RELATIONS 
FOR  CURVILINEAR  THIN  SHELLS 


The  theory  of  thin  curvilinear  shells  depends  on  assumptions  made  in  the  theory  of  elastic 
media.  Specifically,  it  is  assumed  that  the  shell  is  thin  and  that  deformations  are  small.  The 
shell's  mid-surface  is  defined  by  general  curvilinear  orthogonal  coordinates,  0Cj  and  o^,  where  z 
is  the  coordinate  normal  to  the  mid-surface  of  the  element  (figure  4.1).  The  theories  of  thin 
shells  have  been  developed  in  many  books,27'31  and  special  higher-order  theories  are  presented  in 
several  research  reports  and  archival  papers. 


Ui(ocj,  a2,  z)  =  uj(al5  c^)  +  z  p,(al5  a2), 


U2(ap  a2,  z)  =  u2(al5  a2)  +  z  p2(a!,  a2), 


W(ccj,  02,  z)  =  w  (a,,  02)  +  z  w^Oj,  a2)  +  —  w2(Oj,  a2), 

Jm* 


(4.1) 


where  z  is  the  distance  along  the  normal  of  the  mid-surface,  -h/2  <  z  <  h/2. 

The  assumed  displacement  fields  have  in-plane  displacements  Uj  and  u2  and  out-of  plane 
displacement  w,  as  well  as  thickness-shear  deformations  and  p2  and  thickness-stretch 
deformations  w}  and  w2.  In  this  work,  the  thickness-stretch  components  Wj  and  w2  will  be 
neglected  in  applications  for  mid  to  moderately  high  frequencies.  With  the  neglection  of  Wj  and 
w2,  there  are  five  branches  in  the  frequency-wavenumber  spectrum  (two  acoustic  and  three 

optical  branches)  that  extend  the  classical  Love  (thin-shell  membrane-bending)  theory  to  a  much 
higher  frequency  range. 

The  six  elastic  strains  at  any  point  ((Xj,  a2,  z)  are  defined  by  Kraus27  as 


£11  = 


1 


1  +  z/Ri 


(£1  +ZK}), 


e22“I7^i7(52+zK2)l 


E33  =  0, 


Yl2 -T~~7Tr(Yl+ZTl)  +  7'.  At.  (Y2+zt2)» 


1+z/Rj 


1  +  z  /  R2 


Yl3 


Hi 


l-hz/R^ 


LA  9 

Y23  =  I^ 7i? 


(4.2) 
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where  the  quantities  with  a  bar  represent  the  mid-surface  normal  and  shear  strains,  and  the  other 
quantities  reflect  the  linear  deformation  through  the  thickness.  Rj  and  R2  are  the  principal  radii 
of  curvatures  in  the  cq,  a 2  directions,  respectively. 


The  mid-surface  normal  and  shear  strains  are  given  by 


1 

duj 

h  U2 

8hi_  + 

w 

hi 

8a! 

hih2 

8a2 

V 

1 

8u2 

+  _HL 

8h2  + 

w 

h2 

8a2 

hih: 

>  8a  j 

R? 

1 

8u2 

ui 

8hj 

hi 

8a! 

hih2 

8a2  ’ 

1 

8uj 

_ Hi 

8h2 

h2 

8a2 

hih; 

1  8ax ' 

The  independent  shear  deformation  components  of  the  shear  strains  are  given  in  terms  of  (5j  and 
P2  as  follows: 

_  1  8|32  Pi  8hq 
h1h2aa2’ 

T  _  i  ap2  P2  dh2 

h2  8a2  hjh2  8a^ 


where  Tj  and  x2  are  the  twists  associated  with  the  mid-surface,  thus, 

_  1  8w  ui  0 

Mi  =  T-^r__:B-+Pl  ♦ 

hi  c/OC j  Ri 
_  1  dw  u2  D 

^2-m^-r7+P2' 


Let  the  total  twist  x  and  the  total  shear  strain  of  the  mid-surface  y12  be  expressed  by 

x  =  Xj+x2  and  Yi2=Yi+Y2-  (4.6) 
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The  bending  components  of  strain  Kj  and  K2  are  the  changes  of  curvature  of  the  mid-surface  so 
that 


Kj  -  1  dftl  |  P2  dhl  , 
hj  daj  hih2  da2  ’ 

K2  -  1  ^^2  ,  Pi  dh2 
h2  da2  hih2  da! 


(4.7) 
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5.  STRESS-STRAIN  RELATIONS 


The  stress-strain  relationships  for  an  isotropic  elastic  shell  can  be  written  for  the  non¬ 
vanishing  strains  as 

a,.  =  ~  —  2  fell  +ve22]> 

11  l-vz 

a22  =  “ -2^22  +  velll’ 
zz  1  -  vz 

*12  =  GYi2> 

*13=G'Y13=kGY13» 

T23=G'Y23=kGY23>  (51) 


and 


G  = 


E 

2(1 +  v)’ 


K  =  7t2/12, 
G'  =  kG, 


(5.2) 


where  E  is  the  Young’s  modulus,  v  is  the  Poisson's  ratio,  G  is  the  shear  modulus,  G'  is  the 
corrected  shear  modulus  to  compensate  for  the  thickness-shear  deformation,  and  k  is  the 
correction  factor  given  by  Mindlin.  ’ 


15/(16  blank) 


6.  LAGRANGIAN  OF  THE  DYNAMIC  SYSTEM  FOR  A  THIN  SHELL 


The  equations  of  motion  of  the  vibration  of  shells  can  be  derived  from  the  Lagrangian  L 
defined  as 

L  =  -U  +  K  +  X,  (6.1) 

where  U  is  the  strain  energy  of  the  shell,  K  is  the  kinetic  energy  of  the  shell,  and  X  is  the 
potential  for  the  work  done  by  external  surface  distributed  forces. 

The  strain  energy  density  U  per  unit  volume  is  defined  by 

U  =  ^[cyll8ll+a22E22+Yl2t12+Yl3T13+Y23T23]-  (6-2) 

Substituting  the  stress-strain  relations  in  equation  (5.1),  the  strain  energy  density  is  obtained  as 

U=4 — ^“2  [Ell  +  £22  +2ve11e22  +— T~(Yl2  +KYl3  +  kY23)]-  (6-3) 

2  1  -  v  2 

The  kinetic  energy  density  per  unit  volume  is  defined  by 

K  =  ^-ps[Ui +U2+W2],  (6-4) 

where  ps  is  the  density  of  the  shell’s  material,  and  the  dot  indicates  the  partial  derivative  with 
respect  to  time  9/dt. 

A  surface  element  dS  at  any  surface  z  is 

dS  =  (l  +  z/R1)(l  +  z/R2)dS,  (6.5) 

with  the  surface  element  at  mid-surface  dS  defined  as 
dS  =  hih2daidcx2- 

The  total  strain  energy  of  a  shell  of  variable  thickness  h(ax,  a2)  is 
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(6.6) 


h/2  __ 

U=  J  J  J  U(a1,oc25z)dSdz. 

aj  a2  -h  /  2 

The  total  kinetic  energy  of  the  shell  is  thus  given  by 


h/2 

K=  J  J  J  K(a1,a2,z)dSdz.  (6.7) 

al  a2  1 2 

The  potential  for  external  distributed  force  fields  acting  on  the  upper  (h/2)  and  lower  (-h/2) 
surfaces  of  the  shell  is  the  product  of  these  force  fields  and  the  corresponding  displacements  at  the 
upper  and  lower  surfaces,  integrated  over  the  upper  and  lower  surfaces. 

Let  Tu,  >  <lu2 ,  and  q+  be  the  distributed  forces  per  unit  area  in  the  1, 2  and  3  directions, 
applied  on  the  upper  surface  of  the  shell  and  let  q"  ,  q“2 ,  and  q~  be  the  corresponding  forces 
applied  on  the  lower  surface  of  the  shell.  Since  the  surface  element  dS  is  a  function  of  z,  then 


dS+  =  [1  +  — — ][1  +  ]dS , 

2Ri  2R2 

dS'=[l-—  ][1 - — ]dS. 

2Ri  2R2 


(6.8) 


The  total  work  done  by  these  distributed  surface  forces  is  the  product  of  these  forces  and  the 
corresponding  displacements  evaluated  at  the  upper  and  lower  surfaces,  i.e., 


X=  J  J  ((lu1Ul’+<ht2U2  +qwW+)dS+  +  J  J  (qUiUi  +q„2U2 +q~W")dS“  (6.9) 

al  «2  al  a2 

where  the  +/-  superscripts  indicate  the  upper/lower  surfaces  of  the  shell. 

Since  the  dependence  of  the  strains  and  velocities  on  z  is  known  (see  eqs.  (4.1)  and  (4.2)), 
then  one  can  perform  the  integration  over  z  across  the  shell's  thickness.  To  perform  these 
integrations,  expressions  containing  various  terms  in  z  must  be  approximated  as  follows  for 

— « 1,  valid  up  to  the  order  0(h7/R7): 


18 


One  may  interchange  Rj  with  R2,  and  vice  versa.  Substituting  the  expression  for  the  displace¬ 
ment  field  in  equation  (4.1),  then  the  expression  for  the  kinetic  energy  density  in  (6.4)  becomes 

K ==  2  Ps[(ui  +  z$i)2  +  (^2  +  z(32)2  +  w2]. 

Using  equation  (6.10),  the  integral  of  the  kinetic  energy  over  z  only  in  equation  (6.7) 
results  in 
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h/2 


J  KdSdz  =  |p2hj 

-h/2  “ 


1+ 


12RjR2 


2  C 


[uf  +U2  +W2]  +  2(Uj^!  +U2/32)^- 


±+± 

vRi  R2; 


(6.11) 


12 


1+ 


3  h 

20  R2R2 


[A2+A2]|ds. 


The  expression  for  the  total  work  done  by  external  forces  in  equation  (6.9)  can  also  be  simplified 
by  substituting  for  the  displacements  and  the  surface  elements  at  the  upper/lower  surface  so  that 


and 


uf=u!±ipi. 


U|  =  u2±ip2. 


W±  =  w. 


This  results  in  the  following  expression  for  X: 


X"  J  J  (9ujul  +(1u2u2  +<lww  +  mpIPi  +mp2P2)dS, 
«1«2 


where 


qUl  =  (qu,  +qUl)^+ 
qU2  =  (qj,  +qu,/i+- 


,+w 

'  1 

1  > 

J  zl 

V 

J _ 

R2^ 

(qu,  qUl  )> 


qw  —  (qw + qw) 


.2  > 


4RiR2y 


+  — 
2 


hr  i  i  i.  + 


vR1+R2J(qu2  qu2)’ 


.2 


1  +  - 


4RjR2 


+  2^Rj  +  R2  j(q^  qw)’ 


h 

mR  =  — -i 
Pi  2 


(qu,  qUl) 


1+-5f_  +2LT-1  •  1 

4R2R2J  2 


Ri  +Rj(q“i+qui) 


(6.12) 


h 

J2  2 


mp  =-• ! 


(qu,  qu? ) 


i+ 


4RiR2 y  2 


hf  1  1  \  + 

Ri  +Ro  qu2  +qu2) 


VR1 


(6.13) 
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Note  that  if  the  surface  forces  are  applied  equally  at  the  two  surfaces,  i.e.,  if  =  q„2  =  qU[  /  2, 
qu2  =  Qu2  =qu,  /2,  and  q^,  =  q^,  =  qw  /  2,  then  the  force  fields  are  the  resultants  applied  at  the 
central  surface  of  the  shell,  z  =  0,  i.e.: 


Qu]  —  Qu, 


1  +  - 


4RiR2. 


f 


Qu2  “  Qu2 


l1+ 


h2 

4RtR 


A 


2j 


f 


Qw  —  Qw 


1  + 


4Ri  R2 


9 


_  h^ 


mP,  =  T 


f_L+ JL 

^R1  R2 


_  h2 
mP2=<lu2^ 


CL+ n 
R1  R2 


where  q  represent  the  equivalent  force  fields  applied  at  the  central  surface,  z  =  0. 

The  integral  of  the  strain  energy  density  in  equation  (6.3)  can  be  written  as: 

V  =  J  J  J  j(£i  +  zKj)2  A(z)  +  '  +2v(£!  +zK1)(s2  +zK2) 

ccj  a2  -h  /  2  1 


1  —  v 

+ - 

2 


(Yi  +zt1)2A(z)  + 


(Y2  +zt2)2 
A(z) 


+  2(Yi  +zt1)(y2  +zx2) 


+K(pf  A(z)  + 


dzdS, 


(6.14) 


where  A(z)  = 


1+z/R2 
1  +  z/Rj  ' 
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Integrating  the  strain  energy  density  over  z  and  keeping  terms  up  to  0(h3)  in  the  Taylor 
series  expansion,  one  gets 

1  E 


U=5TT^J  JUd5 
«1«2 


where 


U  =  h  jef  +  ej  +  2ve1e2  + ^[(Yi  +  Y2)2  +  k  (jif  +  n2)]j  +  JJ  j(Ki  +  K2  +  ZvKjK^) 


2  2 

Ri~R2~  2£lKl  +  282K2 


Y  i 

_  n 

1  —  V 

-1 

JIri 

R2j 

T 

2 

_  - 

y?  v| 

_2YiTi+2Y2T2 

K1  r2 


vRi  r2; 


+(ti+t2)2+yk 


'"it2  tt2 
i£l_  _  J£2_ 

vri  r2; 


^R1  R2  > 


(6.15) 


Note  that  Yi  +  Y2  =  Yi2-  The  constant  y  in  the  last  term  of  (6.15)  was  introduced  by  Naghdi34  as 
Y  =  3/7  to  further  correct  the  shear  correction  factor,  k. 
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7.  GEOMETRIC  AND  KINEMATICAL  EQUATIONS  FOR  A 
PSS  OF  CONSTANT  THICKNESS 


Consider  a  PSS  of  constant  thickness  whose  mid-surface  is  defined  by  \  -  a  >  1  and  inter- 
focal  distance  d.  A  point  on  the  shell’s  mid-surface  is  given  by  coordinates  0Cj  =  r|  and  ccj  =  <J>. 
The  shell’s  length  2 1  and  maximum  radius  R  are  given  by  (see  figure  7.1) 

Half  length  i  -  da/2, 

Maximum  radius  R  =  dVa2  - 1  /  2, 

Va2  -1 


Ratio  of  diameter/length  = 
Eccentricity  e  =  1/a. 


(7.1) 


The  metric  coefficients  of  a  surface  element  at  the  mid-surface  of  the  shell,  h^  and  h^,  and 
the  surface  element  dA  at  any  point  r|,  <j)  are  given  by 


— 4-  and  h^  =^Va2 -l^l-tl2 , 


and 


dS= 

\2) 


\2 


Va2  -  li/a2  -T)2  dr) d<|) . 


(7.2) 


(7.3) 


The  principal  radii  of  curvatures  at  any  point  (rj,<(>)  are 


Rrj  —  Rl  — 


1  (a2-,,2)3'2 


R4,=R2=r-Va2-lVa2-ri2  . 
Y  2a 


(7.4) 


Letting 

Uj  =  u, 

R1  =  Rtj» 


u2  =  v,  (3,  = 

R2  =  hj  =  h^. 


P2 -  3<t>, 

h2  =  V 
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and 


Ki  =  K,,  = 


j_apri 

h-n 


Ko  =  KA  = 


1  jhj 

h-qh(j)  3t| 


»l  =  \h\ 


M-2  =  f^4> 


1  3w_  _u_  R 

i^3n""5ii  ’ 

1  dw 


Pri» 


h(j)  d<|> 


(7.6) 
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8.  DERIVATION  OF  THE  EQUATIONS  OF  MOTION 


The  equations  of  motion  of  the  shell  are  derived  from  the  Lagrangian.  Let  a  function  L  be 


L  =  L(xj,yk,|^-,  )  =  L(*j,yk,yk,j,yk,ji),  k=  1,2, 3, 4,  5;  jj=  1,2,3,  (8.1) 


where  Xj  are  the  independent  variables,  yk  are  the  dependant  variables,  and  the  comma  in  the 
sub-indices  indicates  partial  differentiation.  The  Lagrangian  L  is  given  as  an  integral,  i.e.. 


F=J  J*  jLdSdt  =  J*  JJ(  -U  +  K  +  X)dSdt. 


«1«2 


al  a2 


The  partial  differential  equations  of  motion  are  then  given  by 


(8.2) 


3l _ a_ 

f  dL  ) 

a  j 

f  dL  ) 

__a_ 

f  dL  ) 

a2 

f  dL  1 

dl 

f  dL  ) 

9yk  dxj 

layk,ij 

3x2 

8x3 

laVk,3j 

9x2 

^yk,n) 

9x2 

^dyk,22  ) 

(8.3) 


9x? 


dL 

a2 

r  dL  > 

a2 

f  dL  ) 

a2 

f  dL  ) 

v^yk,33  j 

j 

3xi 3x2 

l^yk,!2 ; 

9x2  8x3 

^yk,23  ) 

j 

8x3  8xj  1 

,ayk3i> 

=  0, 


where  k  =  1,  2,  3, 4,  5.  Letting 


xi  =r\,  x2  =  4>,  x3  =  t, 

Yl  =  U,  Y2  =  v>  Y3  =  w’  Y4  =  Prp  and  Y5  = 

then  there  are  five  equations  of  motion  on  the  five  independent  kinematical  variables. 

One  can  obtain  the  five  coupled  partial  differential  equations  by  applying  the  calculus  of 
variations,  equation  (8.3),  on  the  Lagrangian  L.  Since  the  Lagrangian  is  a  very  long  expression, 
the  application  of  equation  (8.3)  is  performed  for  each  component  of  L,  i.e.,  K,  X,  and  U. 
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To  simplify  the  algebraic  manipulations,  one  can  multiply  L  by  (l-v2)/(Eh).  To  non- 
dimensionalize  the  independent  variables,  let 


u  =  u /  ^  v  =  \/£  w  -  w  /  £, 


(8.4) 


and 


2  E 


CP  = 


Psd-VZ) 


(8.5) 


where 


da  j^2 

^=-i“  =  half-length  of  the  shell,  and  E  = - tt  is  the  flexure  coefficient. 

2  12 1 


(8.6) 


8.1  EXTERNAL  WORK 

The  external  work  X  has  terms  with  yk  only,  so  that  the  application  of  (8.3)  results  in  the 
following  terms  of  the  type  3X/3yk: 


Eql: 
Eq  2: 
Eq  3: 


1  — 
Eh£ 

1-v2 

EYi£ 

1-v2 

Eh^ 


•qu  h^  h^, 


■qv  h^  h,j,. 


'Tw  ^t|  ^4>  ’ 


(8.7) 


1-v2 

^4:  ij^rmP,Nh+’ 


Eq  5: 


1-v 


Eh^ 


2  N  V 


8.2  KINETIC  ENERGY 


The  inertial  terms  are  determined  by  using  (8.3).  The  kinetic  energy  has  terms  of  the  type 
Yk,t so  that  the  application  of  equations  (8.3)  results  in  differentials  of  the  type  ^ 


at 


,  i.e.. 
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8.3  STRAIN  ENERGY 

The  strain  energy  has  terms  with  partial  differentials  of  up  to  second-order  in  t|  and  t}>. 
Applying  (8.3)  to  the  expressions  for  U  in  (6.15)  would  result  in  five  coupled  partial  differential 
equations.  These  can  be  written  symbolically  as 

Eq  1:  LuuIi  +  Luvv  +  Luvvw  + 

Eq2 

*  vu ^  Lw  V  +  LyWW  +  +  Lyp^P^, 

3.  Lwuii  +  Lwvv  +  Lwww  +  Lwp^p^j  +  Lwp^p(j),  (8.9) 

Eq  4:  Lp^u  +  Lp^v  +  Lp^w  +  Lp^P,  +  Lpj^P* , 

Eq  5:  Lpt„u  +  LpjVv  +  LptWw  +  Lp^  +  Lp^. 

The  partial  differential  operators  (Ly)  will  be  defined  in  the  proceeding  subsections.  To  simplify 
the  expressions,  the  following  variables  will  be  used: 
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A  =  a2  -  1 


B  =  1  -  ri2 


C  =  a2  -  ri2, 


I  a2-n2 


Y  =  3/7,  k  =  7i2/12. 


8.4  EQUATIONS  OF  MOTION  FOR  NONAXISYMMETRIC  VIBRATIONS 


8.4.1  First  Equation 


Luuu  =  DB^ 


,^2x9^  F  92u  2D  2  D  _  2D3 

r)D(l+D  )t~  +  "  ~  2  +  ~~V&Z - Fica2 —  u 

d*1  DB  0(}>2  L  B  C  cj 


+J_a4BD  p3u|1ir  7r>2^U  ,  ta4D3d2u  4  2  D  ^  6D3b1 

1  C*-  B3?  +  ,1(3‘7D  +F— a  ’I  7^-^“  —  4  <810> 


3tiJ  A3  3(f) 


Luvv  =  (1  +  F)-^  +  (v  +  F)^-  +  sj— — j(l  +  F)|4, 
B  a(J)  8ri3(|>  AC2  3(j) 


(8-11) 


Luww  =  ^=^-[(1  +  FK)D2  +  v]g  +  a^-(4A  +  B)w 


(8.12) 


Lup  Pri  =  FKD2P^  +  e(a2^^l  -  4a2^  fa  - 

1  [  C2  3ri2  c3  3t|  AC  dif2 


,  f  a V  ^  a4AB]D  \ 

+[-^—FyK—  M’ 


(8.13) 


L„p,P»=e|-(,+F)a2^ 


(8.14) 


30 


8.4.2  Second  Equation 


L™D = (v + F)l^  ■  i(1 + F)lk E  |"° + F)i?fr 


(8.15) 


ov  9  3v  1  3  V  f  7  D  7  D  9  D~L 

Lvvv  =  FDB— ~--FD(l  +  D2)r|— + - ~-  +  F  a2 - rf - Ka2—  v 

vv  dr]2  dr)  DB8({>2  C  B  AJ 


„  4Bd["_92v  „2xdvl  4 D3  92v  Fa4D|"  2  a2B"|_)  /a  i*\ 

+8 (-Fa4 — s-  B — T  +  ri(3-7D  )—  +a4— T— T - 5 —  r\  +yK— —  vk  (8.16) 

C3  3r)2  3tlJ  A3  8<}>2  C2A  L  A  J 


LvwW_ VAB^(1+FK)+vD^a<t)  +8{(1+FYK)A3/2C2  }’ 


(8.17) 


LvPnPT!-E|(l  +  F)^2A2  M, 


(8.18) 


=  FK^+EL4F^,|i-^^  +  Fa2^ 

v  C3  3ti  AC  0(J)^  Cz  drf 


r7/Ti2a2  a4B^Q 

+  F(-^2-  +  TfK^2)P+ 


(8.19) 


8.4.3  Third  Equation 

Lv™0  = =  -  a  ^r[A(l  +  Fk)  +  vC  ]|5  +  -£L[l  +  (v  +  Fk)D2  -  3Fk  ^  ju 

+e  |  a5  ^§^-(1  +  FkY)^  +  Tl^|.[l  +  3FkyD4(2  -  3D2)]u  j,  (8.20) 

(i + fk + vd2)I + ef (1 + fky)^S}  ’  (8-2i) 
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Lwww  =  FkD 


pd  w  n. Ln2w^  ,  1  d2w 

B— 2""(1  +  D  )f|~r — I — 5 — — ~y 
dr\  dri  D2B  dty2 


a2D 


[l  +  D4  +2vD2]w 


+e 


-Fky 


a4DB 


B^  +  (3-7D2)n^ 


c°  L  3n 


3ri 


a4D3  d2w  i  (,  D3B 

"a3”  a? 


+  Fky — — — o-  +  a 


'J__  1  ■ 

.c3  A3. 

w 

,  (8.22) 


aPv, 


+  —a 


3Vba 


B(Fky  + 1) 


a3VB 

c2Va 


ll[l  +  3FKYD2(l-2D2)PT1]J, 


(8.23) 


,  a  _c  VC  3P*  L  „  ,a3VB3Pj 
wPA  FKaVB  30  +E{[1  +  FKYlI?72-'d^}' 


(8.24) 


8.4.4  Fourth  Equation 
Lp  uu  =  FkD2u  +  e 


f„2  B2  82u  ^  a2AB_c)u  „  a2  d2u  . 

a2ri2  a4AB 

■* 

{  C2  0r|2  C3  AC  9<j>2  + 

— 3 — Fky— t- 
C2  C4 

U  ► 

(8.25) 


(8.26) 


LP„ww  = 


VAB  3w  a3VB 


-Fk- - —  +  8 


a  3ri  C3Va 


r)w 

(1  +  Fky) AB  —  +  nC  (-1  +  3D2  (1  -  2D2 )) w 
dT] 


(8.27) 


LP,p„Pti  ■ 


-FK^P^+e 

a 


Bi%-B(i+D2)|a 

dr|  dr| 


F  a2Pn 
DB  04)2 


+D 


a2B^ 


f  2  2 

n  az  _  „  „  , 

v  B  VC+FKYC2N' 


(8.28) 
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= £  i(v + F)l-lr + (1 + f)“3^- 

W'  I  dria<t»  ^ B  3<J> 


8.4.5  Fifth  Equation 


T  _  L  a2qD3u 

LMu=£{(1+F)^^r 


LP*vv  =  Fkv  +  £ 


Fa2B 

La2v 

A  dv 

a2  92v  F 

2  2  a4B 

c2 

B — T 

L  M2 

-4— ri— 
C  'dr| 

AC  dip2  +  C2 

a  rr  +Ky-^- 

Vc  3w 


a3VB  dwl 


Lr  Ww  =  -Fk — j=— —  +  £(  -(1  +Fky) - Tn 

P*w  aVB  9(j)  1  An 3/2 


AC3/2 


LP«P^  ~8 ) 


(V  +  F)l^_a  +  F)B^Ll’ 

3rj8(J)  B  8(}) 


■FK^4^|3A  +  e 
az 


FDB— ^-Frp(l  +  D2)^-  +  1  -  ^ 
3ri2  9ri  DB  a<J)2 


+  Fa2D 


1  T|  B 

- 4 - KY - 

C  a2B  AC 


8.4.6  Final  Nondimensional  Form  of  the  Equations 
Eq  1:  Luuu  +  Luvv  +  Luww  +  Lyp^P^  +  L^^P^ 


rL  a4  Vac 1  i-v2  Vac 

_c21^  +  8C2  a2  u  +  £^  +  F>  J  ^  ^  ^2  Qu. 


(8.29) 

(8.30) 

(8.31) 

(8.32) 

(8.33) 

(8.34) 

(8.35) 
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Eq  2.  Lvuu  +  Lvvv  +  Lvww+Lvp^P^ +Lvp^ 

4(a  +  4)#f  +  e(  1  +  D^ 


i-v2  .Vac 


Eh  a 


n  v  * 

^  2  ^v’ 


(8.36) 


Eq  3.  Lwuu  +  Lwvv  +  Lwww  +  Lwp  P^+L 


wp, 


.4^a+4)#_i-^.VAc 


Cp  a 


Eh 


?  ‘Iw* 


(8.37) 


Eq  4:  Lp^u  +  Lp^v  +  Lp^w  +  Lp^P^  +  Lp^fy, 


2\s  1  i-v2  Vac. 


(8.38) 


^  5:  LP*uu  +  Lp*vv  +  Lp^w  +  Lp^p^T,  +  Lp^p 

4ej#^+(l  +  D2,vl-i4# 


Cp  l  a- 


eT'^V 


(8.39) 


8.5  EQUATIONS  OF  MOTION  FOR  AXISYMMETRIC  VIBRATIONS 

For  axisymmetric  motion  of  PSS,  a/d(f>  =  0,  and  hence  the  five  coupled  equations  decouple 
into  three  coupled  equations  on  u,  w,  and  Pn  and  two  coupled  equations  on  v  and  (3^.  It  can  be 
seen  that  the  expressions  in  (8.11),  (8.14)  in  the  first  equation,  (8.21)  and  (8.24)  in  the  third 
equation  and  (8.26)  and  (8.29)  in  the  fourth  equation  vanish,  rendering  the  first,  third,  and  fourth 
equations  dependent  on  u,  w,  and  PT1  exclusively.  Furthermore,  for  axisymmetric  motion, 
expressions  in  (8.15),  (8.17),  and  (8.18)  in  the  second  equation,  (8.30),  (8.32),  and  (8.33)  in  the 

fifth  equation  vanish.  This  means  that  the  second  and  fifth  equations  are  exclusively  dependent 
on  v  and 
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9.  REDUCTION  TO  SPHERICAL  SHELLS 


The  equations  of  motion  of  spherical  shells  with  shear  deformations  and  rotatory  inertias 
can  be  recovered  from  section  8  by  taking  the  limits,  d  ->  0,  a  and  -»R  These 

were  also  derived  earlier  by  Wilkinson,35  Prasad,36  and  Wilkinson  and  Kalnins.37 


9.1  FIRST  EQUATION 


Luuu=Vb|^[V§u]  +  F(2-k)u  +  ^|^, 
drf  B  dty 


L^=(I+F)ff+(v+F> 


a2v 

'a^af 


Luww  =  (1  +  v  +  kF)Vb|^-, 


P<t>  =0’ 


and 


Right  side  =^-{(1  +  e)u  +  2e  P^}— 1  Rqu. 


9.2  SECOND  EQUATION 

T  -  ,  a2u  ,,  x,,  ti  9u 

Lvuu  =  (v  +  F)  -  (1  +  F)— — , 

3r|3<f)  B  9<f> 


Lvy7  =  FVB±r[VBv]  +  F(2-K)V  +  I|l. 


Lvww  =  (1  +  v  +  kF) 
lvP,.Pti  =  0* 


1  8w 

Vb 


(9.1) 

(9-2) 

(9.3) 

(9-4) 

(9.5) 

(9.6) 

(9.7) 

(9.8) 


35 


(9.9) 


2  2 
Right  side  =  ^-{(l  +  8)V  +  2s  p^}  — Rq 

C  n  Eh 


(9.10) 


9.3  THIRD  EQUATION 

Lwun  =  -(1  +  V  +  kF)— [Vb  u], 
9ti 

Lwvv  =  -(1  +  v  +  kF)^|, 


.Www  =  Fk— 


B 


dw 

dri 


Fk  d2w 
B  3(()2 


-2(1  + v)w. 


^  wp^  Pri  Ek  ^  [  Vb  ] , 

Lwp,3<(,  =Fk 

2  2 

Right  side  =^2~(l  +  e)W-— ~  V  Rqw. 
cp  Eh 


(9.11) 

(9.12) 

(9.13) 

(9.14) 

(9.15) 

(9.16) 


9.4  FOURTH  EQUATION 

Lpyi  =  FKu, 


LP„vV  =  0, 

L3iiWw  =  -FkVB^, 


lP,P„Ph  =  -FkPii  +  e 


a 

< 

[b^I 

8ri 

.  aTl  . 

l  32P„  n2 
+  F--^-[v  +  4-]pT 


b 


B 


(9.17) 


(9.18) 

(9.19) 


36 


Im<5*=4v+F)S+«+F>stL 


(9.20) 


jk —  r**  T  1 _ 

Right  side  =— j- 8 1^+211 1 — -mp  . 


(9.21) 


9.5  FIFTH  EQUATION 


Lp,»u  =  °. 

LRvv  =  Fkv, 


(9.22) 


T  _  c  1  dw 

lMw  =  "FkVb^’ 


(9.23) 


L^Pti=s(v  +  F)^A-(1  +  F)^., 


(9.24) 


LbbP*=-FkP*  +  e  F^-  B ^  +-^f  +  F  1-J- 

‘w*  *  driL  an  J  b  a<j>2  L  b 


id% 


(9.25) 


R2  To  1-v2 


Right  side  =^-e{P,k+2v} — — ^-mR 
cD  1  J  Eh  K<j> 


(9.26) 


37/(38  blank) 


10.  SUMMARY 


This  report  has  documented  the  derivation  of  the  equations  of  motion  of  PSS  of  constant 
thickness.  The  equations  were  derived  from  the  Lagrangian  of  the  system,  including  the  strain 
and  kinetic  energies  and  the  potential  for  external  mechanical  surface  forces  and  moments.  The 
strain  energy  density  was  derived  for  three  independent  displacements  and  two  changes-of- 
curvature.  The  strain  energy  was  developed  in  a  Taylor  series  of  thickness-to-radius  ratios  of  up 
to  third  order.  Thus,  the  theory  includes  the  effects  of  membrane,  bending,  and  thickness-shear 
deformations.  The  kinetic  energy  density  includes  the  translational  as  well  as  rotational 
components,  so  that  the  equations  of  motion  include  translational  and  rotatory  inertias.  The 
potential  for  external  forces  includes  distributed  surface  forces  and  moments. 

Use  of  Hamilton’s  principle  on  the  Lagrangian  of  the  system  resulted  in  five  coupled  partial 
differential  equations  for  nonaxisymmetric  vibration.  These  equations  are  self-adjoint  and 
positive  definite  due  to  the  use  of  the  Lagrangian  and  the  calculus  of  variations  in  their 
derivation.  The  five  partial  differential  equations  are  coupled  for  non-axisymmetric  motion  but 
reduce  to  two  systems  for  axisymmetric  vibration.  One  system  has  three  coupled  partial 
differential  equations  on  u,  w,  and  (3^,  i.e.,  nontorsional  displacement  fields  and  the  second 
system  has  two  coupled  partial  differential  equations  on  v  and  (3^,  i.e.,  torsional  displacement 
fields. 

The  thin-shell  theory  used  in  this  report  results  in  a  frequency-wavenumber  spectrum  with 
five  branches.  The  accuracy  of  the  thin  shell  theory  is  judged  on  how  accurate  these  branches 
agree  with  the  exact  elasticity  theory.  It  has  been  known  that  the  inclusion  of  the  thickness-shear 
and  rotatory  inertia  effects  makes  the  lowest  flexural  branch  have  the  correct  group  and  phase 
velocities  in  the  high-frequency/high-wavenumber  limit.  The  second  acoustic  branch  represents 
torsional  motion  and  it  also  has  the  correct  phase  and  group  velocities  in  the  high-frequency 
limit.  The  first  optical  branch  represents  longitudinal  motion,  where  the  group  and  phase 
velocities  are  approximately  correct.  The  remaining  two  optical  branches  represent  the 
thickness-shear  modes. 

These  equations  were  reduced  to  those  for  the  nonaxisymmetric  vibration  of  spherical 
shells  and  shown  to  agree  with  those  developed  earlier  by  various  authors. 
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The  five  coupled  partial  differential  equations  have  nonconstant  coefficients  and  do  not 
have  closed  form  eigenfunctions  that  would  satisfy  the  system  of  equations.  In  future  reports, 
approximate  comparison  functions  will  be  used  for  the  five  dependent  variables  u,  v,  w,  |V,  and 
P<t>-  These  would  result  in  coupled  5N  x  5N  algebraic  equations  in  terms  of  N  comparison 
function  sets. 
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APPENDIX 

COMPENDIUM  OF  USEFUL  FORMULAE  FOR 
PROLATE  SPHEROIDAL  SHELLS 


In  this  appendix,  the  following  variables  will  be  used: 
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